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Simulation results are presented on Higgs mass calculations in the spontaneously broken phase of the Higgs 
sector in the minimal Standard Model with a higher derivative regulator. A heavy Higgs particle is found in the 
TeV mass range in the presence of a complex conjugate ghost pair at higher energies. The ghost pair evades easy 
experimental detection. As a finite and unitary theory in the continuum, this model serves as an explicit and 
simple example of a strongly interacting Higgs sector without technicolor. 



1. Introduction 

Following the systematic method that was de- 
veloped earlier [1], we report our simulation re- 
sults on a heavy Higgs particle in the Higgs 
sector of the minimal Standard Model with a 
higher derivative regulator. The quantization 
procedure in Minkowski space-time with indefi- 
nite metric Hilbert space for the 0(4) symmetric 
higher derivative Higgs Lagrangian, 
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(1) 



leads to the euclidean partition function [1,9] 
Z = J X'$exp(— J (I'^xCe), where Ce is ob- 
tained from Eq. (1) by replacing the Minkowski 
operator □ with the euclidean Laplacian. 

The general particle content of the model can 
be exhibited by diagonalizing the free Hamilto- 
nian in terms of creation and annihilation opera- 
tors. The interaction part of the Hamiltonian can 
be expressed in terms of creation and annihilation 
operators as well [9]. 
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2. Unitarity 

It has been pointed out before that field the- 
ories with a complex ghost pair will have a 
unitary S'-matrix in the scattering of ordinary 
particles [8]. The unitarity of the S'-matrix 




Figure 1. The elastic scattering cross section of 
two Goldstone particles is shown in the large N 
approximation as a function of the center of mass 
energy in the isospin 7 = channel . The peak 
corresponds to the Higgs resonance. The scatter- 
ing phase shift of the 7 = 0, J = channel is also 
shown. Note that the ghost effects are hidden and 
the scattering is unitary. 

has been demonstrated within the framework of 
the l/N expansion in the broken phase of our 
higher derivative model [10]. Starting from the 
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Lippmann-Schwinger equation, one can show that 
the elastic scattering amplitude of Goldstone par- 
ticles remains unitary at all energies and the ef- 
fects of the ghost pair are hidden in experimen- 
tally measurable quantities. In Fig.l, a plot of 
the scattering cross section and the correspond- 
ing phase shift in the isospin 7 = and angu- 
lar momentum J = channel of elastic Gold- 
stone scattering is shown. The peak in the cross 
section is the standard Higgs resonance in this 
channel. Surprisingly, as the center of mass en- 
ergy is increased beyond the energy threshold of 
the ghost pair, the cross section remains a very 
smooth function and exhibits no visible signal of 
the ghost pair. This can be understood from the 
fact that the complex ghost mass in our theory 
has a very large imaginary part. If the complex 
mass of the ghost state were tuned very close to 
the real axis, one would expect to see a more visi- 
ble effect on the cross section. Since the S'-matrix 
remains unitary, the Higgs model with a higher 
derivative regulator defines a finite and unitary 
theory in the continuum whose intrinsic cut-off 
scale is not revealed by some wrongly assumed 
pathological properties of the scattering ampli- 
tudes. 

3. Lattice Action and Phase Diagram 

For non-perturbative computer simulations of 
the higher derivative Lagrangian, we introduce a 
hypercubic lattice structure. The lattice spacing 
a defines a new short distance scale with the asso- 
ciated lattice momentum cut-off at A = it /a. We 
will have to work towards the large A/M limit in 
order to eliminate finite lattice effects from the al- 
ready regulated and finite theory. The euclidean 
lattice action that we studied has the form 

Se = [-K <f>aiO + -j^)'t'a + {l-^K)<t)a<t)a 

X 

+A {c^„c^„ - 1)2] , (2) 

where □ is the lattice Laplace operator, with the 
lattice spacing a set to one for convenience. 

The lattice model exhibits two phases, as ex- 
pected. The symmetric phase with full 0(4) sym- 
metry is separated by a second order phase transi- 
tion line from the broken phase which has a resid- 



ual 0(3) symmetry for every fixed value of A in 
the (k, M) plane. 

We developed a Hybrid Monte Carlo algorithm 
with Fourier acceleration which solved the prob- 
lem of critical slowing down very effectively for 
our model [9]. Since the updates are performed 
in momentum space, we can use improved lattice 
propagators without significant cost in computer 
time. This technique greatly extends the param- 
eter range that we can study. Our large N calcu- 
lations indicate that lattice cut-off effects in the 
renormalized coupling constant of the improved 
lattice action are negligible even at M = 2a~^ . 

4. Results and Discussion 

Most of our computer simulations were per- 
formed on lattices which have cylinder geome- 
try in the size range 16^ x 40 to 20^ x 40. On 
each lattice with a given parameter set, 20, 000 to 
60,000 Hybrid Monte Carlo trajectories were ac- 
cumulated. Typical autocorrelation times in our 
simulations are of the order of 10 trajectories. 

A simultaneous fit to rriH and v was performed 
on the data, as described in [9]. The results are 
shown in Table 1 where points marked with * 

Table 1 



K 


M 


A 




niH 


R 


0.056 


0.8 


0.4 


0.057(2) 


0.40(2) 


7.0(4) 


0.056 


0.8 


0.4 


0.045(1) 


0.33(2) 


7.3(5)* 


0.105 


0.8 


0.1 


0.065(2) 


0.31(2) 


4.8(4) 


0.115 


0.8 


0.05 


0.082(1) 


0.24(1) 


2.9(1) 


0.081 


1.0 


0.3 


0.093(1) 


0.42(2) 


4.5(2) 


0.081 


1.0 


0.3 


0.088(1) 


0.38(2) 


4.3(3)* 


0.088 


2.0 


1.0 


0.058(1) 0.351(5) 


6.1(1) 


0.088 


2.0 


1.0 


0.045(1) 


0.29(1) 


6.4(3)* 



correspond to a cylinder size of 20^ x 40. All 
other points correspond to size 16^ x 40. The 
last two points were obtained with the improved 
lattice action while the simple lattice action was 
applied at the other points. 

The radial Higgs excitation in our simulations 
corresponds to a stationary state with real energy. 
Avoided level crossing with the lowest Goldstone 
pair of zero total momentum would only occur in 
larger spatial volumes. However, both the large 
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N calculation and the practical simulation data 
seem to indicate that the ratio niH /f r is not sen- 
sitive to the size of the volume in the range of 
our simulations. A more sophisticated study of 
the Higgs resonance in the finite volume utilizing 
Liischer's prescription is feasible in our model. 
However, preliminary runs show that one needs 
significant computer resources to get reasonably 
accurate results for the scattering phase shift. 

The ghost mass parameter quoted in the table 
is the bare value. It is very hard to determine the 
precise ghost pole location from the Higgs corre- 
lation function due to its mixing with the Higgs 
channel. However, fits to the momentum space 
propagator and analytic calculations in the l/N 
expansion indicate that renormalization effects in 
the ghost mass are not significant. 

The lattice effects of our simulation points were 
also analyzed using the method described in [11]. 
We find that all our simulation points have less 
than one percent Euclidean violation in the prop- 
agator, and they are particularly small for the im- 
proved action. Therefore, the results should rep- 
resent, to a good approximation, the physics de- 
scribed by the continuum higher derivative Higgs 
model. 

It is seen that the niH /f r ratio in our model 
ranges from 3 to 7, and even larger, to be com- 
pared with the earlier ratio of 3, or less [2-6]. 
The Higgs mass range niH > 1 TeV, with a ghost 
location in the multi-TeV region, implies the ex- 
istence of a strongly interacting Higgs sector in 
the model. 

When viewed as a regulator of the conventional 
theory, the effects of the complex ghost parameter 
can be determined in the weak interaction regime. 
The scattering amplitude can be calculated in the 
presence of the ghost pair and compared with 
the corresponding result in conventional pertur- 
bation theory. A comparison can also be made 
within the large N approximation. The large N 
results show that, for a rriH /vr ratio of 3, the ra- 
tio M/niH is of the order of 30 at infinite bare 
coupling. For such a large ghost mass parameter, 
the M dependence in the scattering amplitude is 
practically invisible. Therefore, we have to con- 
clude that the Higgs mass bound at around 700 
GeV in the old calculations was imposed by the 



underlying lattice structure. Heavier Higgs mass 
values can be obtained with the higher derivative 
regulator without any loss of euclidean invariance, 
and the theory can be driven into the strongly in- 
teracting regime. 

This strongly interacting scenario in our model 
is also confirmed by simulations in the symmetric 
phase. We have studied two points in the sym- 
metric phase and measured the renormalized cou- 
pling constant. For the point with a Higgs mass 
of niH = 0.434, at M = 0.8, we get Xr = 4.1(7). 
For the point with niH = 0.352, at M = 2.0, we 
get Xr = 2.2(2). These numbers are to be com- 
pared with Xr 1 for the conventional model at 
comparable Higgs mass values [3,7]. 

The Higher derivative Higgs model is a some- 
what arbitrary extension of the Higgs sector in 
the minimal standard model by introducing a new 
parameter M which represents the threshold of 
new physics. Since the ghost pair evades easy ex- 
perimental detection without violating unitarity, 
or Lorentz invariance, it serves as a toy model 
of the strongly interacting Higgs sector without 
technicolor, a scenario which was excluded in pre- 
vious lattice studies. 
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